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ABSTRACT 

One of the most debated issues in the theoretical modeling of cosmic reioniza- 
tion is the impact of small-mass gravitationally-bound structures. We carry out the 
first numerical investigation of the role of such sterile 'minihaloes', which serve as 
self-shielding screens of ionizing photons. Minihaloes are too small to be properly 
resolved in current large-scale cosmological simulations, and thus we estimate their ef- 
fects using a sub-grid model, considering two cases that bracket their effect within this 
framework. In the 'extreme suppression' case in which minihalo formation ceases once 
a region is partially ionized, their effect on cosmic reionization is modest, reducing 
the volume-averaged ionization fraction by an overall factor of less than 15%. In the 
other extreme, in which minihalo formation is never suppressed, they delay complete 
reionization as much as Az ~ 2, in rough agreement with the results from a previous 
semi-analytical study by the authors. Thus, depending on the details of the minihalo 
formation process, their effect on the overall progress of reionization can range from 
modest to significant, but the minihalo photon consumption is by itself insufficient to 
force an extended reionization epoch. 

Key words: cosmology: theory - galaxies: high-redshift - intergalactic medium - 
radiative transfer 



1 INTRODUCTION 

The study of reionization has recently witnessed an ex- 
plosion of observational results. Since the advent of the 
new century, mor e than 10 quasars with z > 5.7 have 
been discovered JWhite et all 120031: iFanl 12004) . indicat- 
ing a strong increment of the neutral Lyman-Q optical 
depth ijOunn fc Peterson! llQril GP) with increasing red- 
shift, often interpreted as the end of reionization. At 
the same time observations of Cosmic Microwave Back- 
ground (CMB) anisotropics have constrained the total op- 
tical depth of CMB photons to Thomson scattering in the 
range r e = 0.17 ± 0.04, where the quoted un certainty de- 
pend s on the analysis t echnique employed l|Kogut et alJ 
120031 ISnergel et al"H2003h . Independent of the details, this 
optical depth can be converted into a reionization redshift 
that is ^ 10. These tests,together with other probes of the 
high-z universe as, such as Lya emission and Gamma Ray 
Bursts (see Ciardi & Ferrara 2005 and references therein), 
provide an invaluable set of observational information. But 
only measurements of the 21 cm line emission from the neu- 
tral intergalactic medium (IGM) and collapsed haloes will be 



able to map the temporal and spatial evolution of the process 



itself (e. o.LMadau. Meiksin fc Reed 


1997tlShaver et alJll999l; 


Tozzi et alJ l200d: llliev et alJ 12002 


. 120031 ICiardi & Madaul 


2003:lFurlanetto. Sokasian & Hcrnouist 2004). 



Numerical simulations and semi-analytical studies have 
strug gled to match this early reionization (for a review 
see iBarkana fc Loebl 1200 ll : ICiardi fc Ferraral I2005T) , which 
requires an enhanced ionizing photon emission at high 
redshift, while delaying complete overlap until z ~ 6. 
Some models have even included scenarios such as dou- 
ble reionization through very massive, metal-free stars 
at high re dshift , and more normal stars at lower red- 
shift (e.g. |Cen| l2003t IWvithe fc Loebl l2003h . or early 
partial reionization due to a decaying particle or X- 
ray photons, fol lowed by l ater full-reionization by nor- 
mal stars (e.g. loiil l200ll: IChen fc Kamionkowski liool 
lHansen fc Haimanl 12004 iRicotti fc Qstrikerl 12004) . More 
conservatively. lciaML^^rrai^^^WiruTl2003l CFW here- 
after) have shown that normal stars with a slightly top- 
heavy Initial Mass Function (IMF) could suffice to pro- 
duce an early reionization epoch, in agreement with the 
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observations of the Wi lkinson Microwave Anisotropy Probe 
(WMAP) satellite fe.o. lKogut et al.ll2003l1 . 

The basic conclusion one can draw from these mod- 
els is that boosting the ionizing photon production to the 
level required for early reionization is not a problem. Nev- 
ertheless, some tension might arise when the models are 
combined with the results from the Gunn-Peterson effect, 
as the reionization history must be tuned to (i) produce a 
large optical depth (to comply with the WMAP constraint), 
and (ii) to complete reionization at a relatively low red- 
shift ~ 6 (as suggested by the GP data). Note, however, 
that this tension might be simply an artifact of the limited 
statis tics of QSO absorption lin es at these high redshifts 
(e.Q. IChoudhurv fc Ferrarall2005l) . A scenario satisfying the 
above requirements could be one in which cosmic reioniza- 
tion starts at very high redshift and then goes through a 
long phase of partial ionization lasting approximately up to 
z ~ 6. 

Such a long-lasting phase might be caused by the 
presence of small-scale, gravitatio nally-bound structures 
lllliev. Scannapieco fc Shapiroll20o3 . hereafter ISS). In hier- 
archical theories like the cold dark matter (CDM) model, 
the smallest structures are the first to collapse and viri- 
alize, and in order to form stars they must radiate their 
virial energy. However, in a purely atomic gas of primordial 
composition, radiative cooling is ineffective below 10 4 K. At 
these temperatures molecular hydrogen is the main coolant. 
But H2 is easily dissociated by UV photons in the Lyman- 
Werner bands between 11.2 and 1 3.6 eV, which are copi- 
ously produced by the firs t stars f e.g. lHaiman. Rees. fc Loebl 
Il997l: ICiardi et al. II2OO0I) . On the other hand, an early X- 
ray background could provide enough free electrons to pro- 
mote H2 formation. Although the relative strength of the 
above feedback effects is controversial, the positive feedback 
from X-rays would not be able to b alance UV photodisso- 
ciation in the smallest obje cts (e.g. lHaiman. Abel fc Reesl 
120001: iGlover fc Brand! liboa^ . Thus, we can expect a pop- 
ulation of sterile "minihaloes" (MHs), which serve as self- 
shielding screens of ionizing photons and increase the local 
effective recombination rate of the gas. This prolongs the 
partial reionization phase and delays the end of the reion- 
ization process. 

When an intergalactic I-front encounters a minihalo 
during cosmic reionization, the minihalo traps the I-front, 
converting it from a supersonic, weak R-type front to a 
subsonic, D-type front, which photoevaporates the mini- 
halo gas. The first detailed studies of this interaction, using 
numerical gas dyna mics simulations with radiat ive trans- 
fer, were r eported in[Sh^iro^iev fcR^tgal i2004L hereafter, 
SIR) and llliev. Shapiro, fc Rasa I fcOOa. hereafter. ISR). 
Since the gas in MHs is much denser than in the sur- 
rounding IGM, the recombination rate during this process 
is similarly increased. Furthermore, MHs were so abun- 
dant during reionization that a sight-line between any 
two given ionizing sources would have passed through 
many of them. This means that MHs have the poten- 
tial to trap intergalactic I-fronts before they overlapped, 
and suggests that they may have drastically increased 



the global consumption of ionizing 


; photons during reion- 


ization (Haiman, Abel & Madau 


2000 


; Barkana & Locb 


pool IShapird l200ll 


IShapiro. Iliev. Raea & Martell 120031 


Shapiro. Iliev & Raeall 


12004 llliev. Shapiro. & Raea 1120051). 



While approximate schemes to account for the effect of 
MHs on photon consumpti on during reionization had been 
previously proposed (e.g. lHaiman. Abel fc Madaul l2000l : 
iBarkana fc Loebll2002h . they did not have the advantage of 
incorporating the results of the detailed numerical simula- 
tions of ionization front-MH interactions reported in SIR 
and ISR. These last studies not only highlighted the impor- 
tance of accounting for the gradual "peeling away" of neutral 
MH gas by ionization fronts, but also provided detailed fits 
to the total number of ionizing photons absorbed as a func- 
tion of MH mass, source flux level, and redshift. These fits 
then serve as a sound founda tion for more detailed analyti- 
cal treatments (such as llliev. Scannapieco fc Shapiroll2005h 
and numerical investigations, that is the subject taken up 
here. 

In the semi-analytical treatment in ISS, the rate of ex- 
pansion of the I-fronts around individual source halos is 
modified to include the effect of MHs. This was done by gen- 
eralizing the I-front continuity jump condition from Shapiro 
and Giroux (1987) to include the consumption of extra ioniz- 
ing photons by MHs, as computed in SIR and ISR. The time 
evolution of the total ionized volume fraction of the universe 
was then calculated by summing these time- varying ionized 
volumes created by each source, over the statistical distri- 
bution of source halos in the universe. It was found that 
MHs increased the number of ionizing photons required to 
complete reionization by as much as a factor of 2, delaying 
completion of reionization by a redshift interval Az ~ 2. 
The fact that the mass fraction collapsed into MHs grows 
over time, suggested that the extra photon consumption by 
MHs would be more efficient at slowing the advance of I- 
fronts at late times. Instead, we found that the spatial clus- 
tering around each ionizing source kept the density of MHs 
fairly constant, where they were actually encountered by the 
I-fronts. As a result, MHs decreased the availability of pho- 
tons to ionize the IGM at all epochs, delaying reionization, 
without extending its duration sufficiently to reconcile the 
GP and WMAP constraints. However, a more significant ef- 
fect in that direction was found when the rising value of the 
small-scale clumping factor of the IGM was also taken into 
account . 

Here we focus on a complementary approach to this 
problem, in which the extra photon consumption by MHs is 
incorporated as a "subgrid" correction in a numerical sim- 
ulation of reionization. We use N-body results to identify 
the ionizing sources and the evolving IGM gas density field, 
a semi-analytical approach to identify the local density of 
MHs, and numerical radiative transfer techniques to prop- 
agate the I-fronts in three dimensions. The advantage of 
the numerical treatment is that it is liberated from the as- 
sumptions of spherically-averaged I-front propagation and a 
reliance on analytical models for clustering and bias. The 
numerical treatment, for example, is better equipped to ac- 
count for the effects of clustered sources, which can result in 
H II regions powered by multiple sources, either simultane- 
ously or sequentially. On the other hand, the semi-analytical 
approach has the advantage that it is not limited by finite 
numerical resolution and dynamic range, while these issues 
are inherent in the numerical treatment. Our hope is that, 
by offering these two distinct approaches, we can gain confi- 
dence in their answers, to the extent that they agree, remain 
cautious, to the extent that they do not, and gain more in- 
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sight into the important effects than either one affords, on 
its own. 

The main goal of this paper is to assess if the photon 
consumption provided by MHs can result in a slow and grad- 
ual reionization, thus offering a natural explanation for the 
high value of r e and, at the same time, account for the late 
reionization epoch hinted by the GP opacity. The paper is 
organized as it follows. In § [5] and [3] we describe the sim- 
ulations of cosmic reionization in the absence of MHs and 
the physics of MHs photoevaporation, respectively. In §0]we 
discuss the implementation of photoevaporation physics in 
the simulation of reionization. In §|S]we present our results 
and in §|5]we summarize our conclusions. 



2 SIMULATIONS OF COSMIC REIONIZATION 

In this study we re compute the simulations of cosm ic reion- 
ization described in ICiardi. Stoehr fc Whitel <2003t hereafter 
CSW) and CFW, now including the effect of unresolved 
MHs. In this Section we briefly summarize the main fea- 
tures of the simulations. We refer the reader to the above 
papers for more details. 

The simulations, performed with t he N-body code 
GADGET dSpringel. Yoshida fc Whitel 1200 J) , follow the evo- 
lution of an "average" regio n of the Universe 1 . The 
"re-s imulation" technique (e.g. iTormen. Bouchet fc White! 
Il997l) has been used to follow at higher resolution the dark 
matter distribution within an approximately spherical re- 
gion of diameter ~ 50/i -1 Mpc included in a much larger 
volum e (479/i -1 Mpc on a side: lYoshida. Sheth fc Diaferiol 
l200fT) . The position and mass of dark matter haloes 
were determined with a friends-of-friends algorithm, while 
gravitationally-bound s ubstructures within t he haloes were 
identified with SUBFIND iSpringel et al.ll200fr) and were used 
to build the merging tree for haloes and subhaloes. The 
smallest resolved haloes (which start to form at z ~ f 9) have 
masses of M ~ few x 10 9 Mq, and the gala xy population was 
model ed via the semi-analytic technique of lKauffmann et alJ 
(1999). A catalogue of galaxies for each of the simulation 
outputs was obtained, containing for each galaxy, among 
oth er quantities, its position, mass and star-formation rate 
(see lStoehrll2003h . 

Within the high-resolution spherical sub-region, a cube 
of comoving side L = 20/i" 1 Mpc was extracted to study the 
details of the reionization pr ocess, using the Monte Carlo ra- 
diative transfer co de CRASH dMaselli. Ferrara fc CiardilEooll 
ICiardi et al.ll20o"ll . CFMR hereafter) to model the propaga- 
tion of ionizing photons into the IGM. 

The only inputs necessary for the calculation are the 
gas density field and ionization state, as well as the source 
position and emission properties, which are provided at each 



1 Throughout this study we assume a ACDM cosmology with 
(n m ,tt A ,tt b ,h,a$,n) = (0.3,0.7,0.04,0.7,0.9, 1), where fi m , fi A , 
and Hf, are the total matter, vacuum, and baryonic densities in 
units of the critical density, h is the Hubble constant in units of 
100 km s -1 Mpc -1 , erg is the standard deviation of linear density 
fluctuations on the 8fe — x Mpc scale at p resent, and n is the index 
of the primordial power spectrum (e.g. [Spergel et alj| 2003[) with 
the t ransfer function taken from Ef stathiou. Bond. &: White I 
(1992). 



output of the galaxy formation simulation described above. 
In particular, the dark matter density distribution is tab- 
ulated on a m esh using a Triangular Shap ed Cloud (TSC) 
interpolation IIHocknev fc Eastwoodlll98H) . Assuming that 
the gas distribution follows that of the dark matter, the gas 
density in each cell grid is set requiring that flj, — 0.04. A 
number N c — 128 3 of cells have been used. To reduce the 
computational cost of the radiative transfer calculation, all 
the sources inside a grid cell are grouped into a single source 
placed at the cell center. Mass and luminosity conservation 
are assured. 



3 MINIHALO PHOTOEVAPORATION 

In the simulations described in the previous Section, the 
smallest resolution element is a grid cell, 156 h~ kpc on a 
side. Minihaloes are thus well below the resolution of the 
N-body simulation, and there are up to thousands of them 
inside each cell. The aim of this paper then is to model the 
presence and photoevaporation of MHs as sub-grid physics. 

According to ISR, the total number of ionizing photons 
absorbed per minihalo atom during photoevaporation can be 
expressed as a function of the minihalo mass, M, redshift, 
z, and level of external ionizing photon flux, F, as follows 



£mh(M,z,F ) = l + cj> 1 ( M -r)Mz)M F o), 
where M 7 = M/(W 7 M Q ), and 

F 



{10 56 s-VK^iMpc) 2 ]}' 



(1) 



(2) 



For sources with a 5 x 10 4 K blackbody spectrum, typi- 
cal for O-stars, 0i(M) = 4.4 (M 7 °- 334+0 - d23 lgM7 ), <j> 2 (z) = 
(l + z)/10 and fo{F ) = F °i99-o-042i g F _ By combining 

equation Q with the average number of minihaloes in a 
given volume, we can compute statistically the average ion- 
izing photon consumption by minihaloes per total atom in 
this volume as 



Pm 



Aima> 



dM^^lMaM, Z ,F ), 



(3) 



where p m is the mean cosmological mass density and 
dn(M,z,S)/dM is the number density of minihaloes in a 
region with an overdensity S = p m /p m — 1. Here the up- 
per mass limit is set by the requirement that MHs are 
not able to cool atomically, that is IW,max = 10 4 K or 
M max = 2.8 x 10 9 (1 + z)~ 3/2 M© in our cosmology. Simi- 
larly, the lower mass limit, M m i n , is taken to be the (instan- 
taneous) Jeans mass in the cold neutral medium. Note that 
we neglect the finite time delay for the gas to r espond hy- 
drodynamically to cooling JCnedin fc Huilll996t) . Account- 
ing for this effect would raise the effective "filtering scale" 
in neutral regions, decreasing the impact of minihalos some- 
what . 

Adopting an extended PS approach llLacev fc Cold 

1993), we approximate the biased number density of mini- 
haloes, dn(M, z, 8)/dM , as 



dn 
~dM 



prn 

~M 



da 2 (M) 



dM 



f[l.68D- 1 (z)-5 L (8),a 2 (M)], (4) 



where D(z) is the linear growth factor at a redshift z, a 2 (M) 
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is the variance of linear fluctuations within a sphere contain- 
ing a mass M, and 



Ah 



f(SL,a 2 



L / UL 



(5) 



Finally, Sl(8) is the linear overdensity corresponding to a 
nonlinear overdensity of S. If 5 > these quantities can be 
related by the standard top-hat collapse model in terms of 
a "collapse parameter" 9 as: 



. 9 (0- sine) 2 1 
o = ttttt - 1, 



2 (1 - cos 0) 3 



and 



5l = 



o /Qs 2/3 

(I) (0-sm0) 



2/3 



(6) 



(7) 



In underdense regions, S and <5l can b e related by the so- 
lution given by iHeathl il977t) (see also iFriedmann fc Piranl 



20 

8 = V (Sl) 
where 



1, 



(8) 



= —j + exp[-0.01 (4.73<5| + 0.83<5| + 0.10(5!)], (9) 

is the ratio of the comoving size of a perturbation to its ini- 
tial comoving size. Note that eq. @J is a Lagrangian number 
density, which is normalized by the mean density p m in eq. 

rather than an Eulerian density, which would have been 
normalized by (S + l)p m instead. 

Using eqs. we can then re-write equation (|3J as 

i(z,S,Fb) = /co11,mh(z, 5) + ^> 3 (F )I(z, S), (10) 
where 



fc 



,^dn(M,z,5) 
dM — , ' M, 
dM 



P m Jm„ llll 

is the collapsed mass fraction in minihaloes and 



(11) 



I(z,S) = — 

P™ J M„ 



dM M^l MMM)Mz) , (12) 



such that <j>a(Fa)I(z,fi) is the extra number of ionizing pho- 
tons absorbed per minihalo hydrogen atom. 



,MH, 



(13) 



where p m is the gas density in the cell and Al is its 
linear physical dimension. For our reference simulation 
/co11,mh(2, 6) is calculated at (z cr ,8 cr ), corresponding to the 
first photon package crossing of the cell and never increased 
thereafter, i.e. once a cell has been crossed by photons and 
fully or partially ionized, no new MHs are allowed to subse- 
quently form in the cell. In almost all cases, complete ion- 
ization of a given cell occurs within 10 6 years of the initial 
illumination. We call this reference case 'extreme suppres- 
sion' and it places a lower limit to the effect of minihaloes 
on reionization. 

On the other extreme, we consider the case in which we 
allow MHs formation in partially ionized and recombined 
regions. In this case, / co ii,mh is updated at each time step 
and it never becomes zero (we call this the 'no suppression' 
case). This is meant to put an upper limit on the effect of 
MHs, but the formation of sterile MHs in highly ionized re- 
gions is not supported by any strong physical motivation. 
In any case, although the impact of feedback on the forma- 
tion of primordial, small -mass objects has been investigated 
by several authors (e.g. lGnedinll200ol : lOh fc Haimar]l2003l : 
lO'Shea et jD EflOj) 

a consensus has not yet been reached. 
In fact, an early X-ray background would heat the IGM be- 
fore ionization from softer photons occurs and may provide 
an entropy floor even stricter than our "no suppression" 
case. As the presence and effectiveness of t his mechanism 
rema i ns highly uncertain how ever (see e.g. lOh fc Haimanl 
120031: iKuhlen fc Madaul I2005T) . we can cautiously consider 
our calculations as providing reasonable estimates of as up- 
per and lower limits on the impact of MHs on cosmic reion- 
ization. 

To estimate the fraction of ionizing photons absorbed 
by MHs, we evaluate the optical depth of a given cell due to 
minihaloes alone, tmh, as: 



Tmh = 0.56<t(i/)7ih/coU,mhAL 



(14) 



Here cr(v) is the photoionization cross section at a frequency 
v and the factor / = 0.56 accounts for the average photon 
path length through a cell (CFMR). The total optical depth 
of a cell includes the contribution of both the diffuse gas, 
r dift , and MHs, r = r dift + t M h, with: 



Tdiff = 0.56<7(i/)nm(l — fa 



i)Al. 



(15) 



4 REIONIZATION SIMULATIONS WITH 
MINIHALOES 

In this Section we discuss the implementation used to in- 
clude the physics of the MHs (§ into the simulation of 
cosmic reionization (§EU- In general, the galaxy formation 
simulations and the radiative transfer calculations are per- 
formed as described in § 3 of CSW, and only some of the 
quantities are modified to include the mean contribution of 
MHs, as follows. 

The gas density field provided by the galaxy formation 
simulations is split into a diffuse component (the IGM) and 
MHs, according to the fraction of gas collapsed into MHs, 
/coii.MH feauation llll . That is, if tih is the hydrogen number 
density in a given cell, only a fraction nii(l — /coII.mh) is 
assigned to the diffuse gas component. From / ca ii,MH we can 
derive the total mass of MHs in the cell as: 



In each cell along the path of the photon packet, the value 
of the hydrogen ionization fraction is updated according to 
eq. (10) of CFMR. In this case, though, the number of pho- 
tons deposited in the diffuse component of the cell is re- 
placed by A r 7 exp(— tmh), to take into account the photons 
absorbed by the MHs. Here JV 7 is the number of photons 
contained in the monochromatic packet that illuminates the 
cell. Any photons in the package not absorbed by either the 
diffuse gas or the MHs in the cell are passed onto the next 
cell. 

To determine when the MHs in a given cell should be 
considered completely photoevaporated in the 'extreme sup- 
pression' case, we proceed as follows. If (A^ 7i mh)i is the num- 
ber of photons absorbed by MHs in a given cell at time step 
i of the simulation, the fraction of MHs mass that these 
photons photoevaporate is: 
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MH,cvapJz — 



(JV 7 ,mh)i 



(16) 



where £ is given by equation 1101 . The mass evaporated 
from the minihaloes, Mmh^mh.cvip, is transferred to the 
IGM and / co ll,MH is modified accordingly. Note that, as we 
assume that once the cell has been illuminated no new MHs 
form, at each step only the flux dependent coefficient in 
equation IIOL (j>s(Fo), is updated according to the value of 
the current flux reaching the cell, F, while the rest is kept 
constant at (z cr ,5 cr ). F is given by: 



F = 



N-. 



(17) 



where At is the time elapsed since a photon packet has 
gone through the cell and d is the distance from the source. 
Once the sum y^.(J-MH,evap)i becomes equal to unity, we set 
/co11,mh = and consider the minihaloes gone thereafter. 
Only the number of photons necessary to complete photoe- 
vaporation are absorbed during the final step in which this 
condition is met. 

In the 'no suppression' case instead, formation of mini- 
haloes can occur at any time. For this reason, all the quan- 
tities in equation 1101 are always updated. In particular, at 
each time step, we first calculate how much new mass has 
gone in MHs and update / co ii,mh(z, 8) accordingly. We then 
subtract the evaporated mass from the MHs and transfer 
it to the IGM. In this case, when the sum ^.(^MH^vapJi 
becomes equal to unity, new MHs are still allowed to form. 

Finally, a slightly different approach is used for MHs 
inside a computational cell containing an ionizing source (see 
Appendix) . 
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Figure 1. Evolution of the volume-averaged ionization fraction 
for three different simulation runs (see text for details) : run with 
no MHs (triangles), run with MHs and 'extreme suppression' of 
their formation (open circles), and run with MHs and 'no suppres- 
sion' of their formation (filled circles). Results are for an 'early' 
reionization case. 



5 RESULTS 

5.1 Reionization histories 

We quantify the impact of MHs on the progress of cosmic 
reionization by computing the volume-averaged ionization 
fraction, function of redshift, as shown in Figure 

In our reference run (open circles) we adopt an emission 
spectrum typical of metal-free stars, a mildly top-heavy Lar- 
son IMF with characteristic mass of 5 Mq and an effective 
escape fraction of ionizing photons f esc — 20%. 

From a comparison with the analogous run in the ab- 
sence of MHs (triangles; L20, 'early' reionization case in 
CFW), it is clear that the presence of MHs reduces the 
average ionization fraction most at early times, where the 
difference is m 15%, while at later times the difference is 
10%. This reflects the fact that less than 20% of the emitted 
photons in our model are absorbed by MHs at high redshifts, 
and this fraction decreases with decreasing redshift. Of the 
above percentage only ,$ 5% is absorbed by MHs in the 
source cells. This is mainly because there are only ~ 150 
(4000) source cells at z ~ 19 (10), while there are about 
2 x 10 6 cells in the simulation. In both cases reionization is 
complete by z ion ~ 13. 

In Figure |21 we show the evolution of the MH collapsed 
fraction (ignoring photoevaporation, dotted line, open cir- 
cles) and the fraction of photoevaporated MH mass (dashed 
line, open circles). The MHs collapsed fraction increases with 
decreasing redshift, until it reaches a plateau, corresponding 



to the epoch at which all the cells in the simulation have been 
crossed by photons and no more MH formation is allowed. 
Note that the total mass in the simulation box should be 
constant, but small fluctuations are possible as it is cut from 
a spherical region that is embedded in a much larger volume. 
The fraction of photoevaporated mass of MHs increases with 
time, until a value of the order of unity is reached at z ~ 14. 
Note that reionization of the IGM and complete MHs pho- 
toevaporation do not necessarily coincide, as there could be 
cells in which the ionization fraction of the diffuse gas is ~ 1, 
but some MHs still survive photoevaporation, or vice- versa. 

The 'no suppression' case (filled circles in Fig. shows 
a much stronger effect due to MHs, which absorb 80-90% 
of the ionizing photons. As a result, reionization is delayed 
by Az ~ 2. These results are similar to those obtained in 
ISS, which was largely focused on the impact of MH around 
individual sources, which were later summed together to give 
a rough estimate of the total MH impact on reionization. 
This means that each source in ISS saw a full compliment 
of MH, a picture very much like the 'no suppression' case 
studied here. 

In Figure [5] the fraction of mass in MHs (dotted line, 
filled circles) follows a trend similar to the 'extreme suppres- 
sion case', but now the abundance of MHs increases at late 
times, and at early times the fractional abundance of MHs 
is slightly higher. At all redshifts, the fraction of mass in 
MHs that is photoevaporated (dashed line, filled circles) is 
much smaller and increases more slowly than the 'extreme 
suppression' case. This is because MHs are continuously 
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Figure 2. Redshift evolution of the fraction of total mass that is 
in MHs (excluding photoevaporation, dotted lines) and the frac- 
tion of mass in MHs which is photoevaporated (dashed lines). 
Filled (open) circles indicate a case of 'no suppression' ('extreme 
suppression') of MHs formation. Results are for an 'early' reion- 
ization case. 



re-formed, while in the previous case no re-formation fol- 
lows the almost instantaneous photoevaporation. As a con- 
sequence of re-formation, MHs survive also after complete 
reionization of the diffuse gas. If this is the case, their pres- 
ence would still be visible in terms of, e.g., a contribution 
to the optical depth to ionizing photons and to the 21 cm 
emission line. 

Finally, we have also explored the 'late' (S5) reioniza- 
tion case of CFW, in which a Salpeter IMF and f e3C — 5% 
are adopted. Carrying out an 'extreme suppression' simu- 
lation, we find that the effect of MHs is qualitatively and 
quantitatively similar to our 'early' reionization model in 
which no MH reformation was allowed. Less than 20% of 
photons are absorbed by MHs at any redshift. This may re- 
flect the fact that our MH collapsed fraction did not increase 
much with decreasing redshift above the level obtained at 
the higher redshifts of our early reionization case plotted in 
Figure Here we note that the MH collapsed fraction is 
somewhat lower for this study than it was in ISS, due to our 
choice of the unknown small-scale power spectrum and, to a 
lesser extent, our coarse-grain density field. We discuss their 
impact on our results in Sec. 5.3. 



5.2 Thomson Scattering Optical Depth 

From the above reionization histories it is possible to derive 
the expected optical depth to electron scattering, r e , as: 



where ar = 6.65 x 10~ 25 cm 2 is the Thomson cross section 
and n e (z') is the mean electron number density at z . Prior 
to complete reionization, n e is obtained from the simula- 
tions at each redshift as n e — y^ Nc XjUijl — fc6U.,MH)/N c , 
where the sum is performed over all the cells and Xi is the 
ionization fraction in cell i. Once reionization is completed, 
we simply assume complete hydrogen and He I ionization 
throughout the box. We also assume complete He II reion- 
ization after z = 3. We find that the Thomson scattering 
optical depth is r e = 0.16 for the 'early' reionization case 
without MH reformation, r e = 0.12 for the 'early' case in 
which re-formation is allowed, and r e — 0.10 for the 'late' 
case without re-formation. 



5.3 Model Uncertainties 

In this Section we discuss the approximations and assump- 
tions we adopted to model the effect of MHs on the progress 
of reionization. 

The MH photoevaporation simulations by ISR on which 
we based our model span a range of fluxes 0.01 < Fo < 1000. 
Very close to the ionizing sources, however, the normalized 
photon flux Fo in equation 11711 can exceed 1000. In the 
framework of our galaxy model this occurs rarely, in less 
than 2% of the cases, and thus it affects only a small frac- 
tion of the volume of the simulation. In such cases we set 
Fo = 1000 in equation fT| . which is a reasonable approx- 
imation since at high fluxes the MHs photon consumption 
dependence on the external flux becomes relatively weak. 
However, under different assumptions about the sources life- 
times, e.g. short-lived sources, the photon consumption by 
minihaloes would be higher. In a higher percentage of cases 
(~ 10%), Fo becomes smaller than 0.01. When the flux be- 
comes so low, relatively few additional photons per atom 
are needed to evaporate the MHs and thus we assume that 
just one photon per atom is absorbed if Fo < 0.01. Both of 
these approximations are small, but conservative in terms of 
estimating the photon consumption by MHs. 

For each output of the N-body simulations, the dark 
matter density distribution (and consequently the gas dis- 
tribution, which is assumed to follow that of the dark mat- 
ter) is tabulated on a mes h using a TSC interpolation 
jHocknev fc Eastwood! Il98ll) . The TSC interpolation has 
the advantage of producing smoothed density values at the 
positions of the grid cells, as required for the radiative trans- 
fer computation. However it is not adaptive and density 
variations on sub-grid scales are thus smoothed out, reduc- 
ing the contrast of the dark matter density distribution. 
In CSW this effect is discussed in terms of a comparison 
between the TSC interpolation and a 32-particle SPH ker- 
nel, which should better capture the density contrasts. They 
found that the overall agreement for low and intermediate 
density is reasonable, but the smoothing of the TSC scheme 
is clearly visible at high densities, especially at the lower 
redshifts (Fig. 9 of CSW). This can lead to an underesti- 



mate of the value of 8 at 



12 and dark matter densities 



p £ 10 (3 4) M Q /ft(kpc/ft) 3 , depending on z. To estimate 
the error introduced by the TSC scheme, we computed the 
ratio of the average number of extra photons absorbed per 
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unit hydrogen atom at z = 10 in both the TSC and SPH 
scheme. This is 

Er^(Fo T f c )p^(io)/(io,^ Tsc ) fl9) 

E^3« H )p m (10)7(10,C H )' 

In general, <)>3(FqJ c ) and ^(-F^J 1 *) need not be equal and 
will depend in detail on the formation and propagation of 
the ionization fronts in the simulations. However, as a rough 
estimate we simply factor out these terms, leaving a ratio of 
1.030, or an error of ~ 3% due to smoothing. 

Due to our coarse-grained density field the total mini- 
halo collapsed fraction in our computational volume ob- 
tained using our model (Fig. |5J is slightly lower than the 
global collapsed fraction obtained directly from analytical 
estimates (e.g. using PS or ST mass functions), particularly 
at later times when density fluctuations become more non- 
linear. Higher-resolution N-body simulations would be re- 
quired to better calibrate the mass function-local overden- 
sity relationship better, but this goes beyond the scope of 
this paper. Note that these effects are most severe in cells 
containing ionizing sources, in which internal density con- 
trasts are higher. 

The semi-analytical treatment by ISS found that the 
importance of MHs as consumers of ionizing photons was 
higher if each source emitted its total lifetime supply of ion- 
izing photons in a short burst at higher luminosity than if it 
emitted the same number of photons continuously over time 
at lower luminosity. In the Monte-Carlo reionization simu- 
lations reported here, each source is assumed to release its 
supply of ionizing photons spread out over the time between 
time-slices of the density field which were provided by the 
galaxy formation simulations. This results in longer source 
lifetimes and smaller luminosities than the fiducial case con- 
sidered by ISS, involving short-lived sources. It is possible, 
therefore, that the MH correction effect would have been 
higher if we had, instead, assumed shorter source lifetimes. 

As there are currently no direct constraints on the 
power spectrum of density fluctuations at minihalo scales, 
our estimates require a large, uncertain extrapolation. Al- 
though the transfer function we adopted is consistent with 
the simulations, assuming a different transfer fun ction, 
such as the often-used one bv lEisenstein fc Hul Jl999t) . can 
change the minihalo numbers significantly. For this particu- 
lar application, the coll apsed fraction in minihaloes for the 
lEisenstein fc Hul (|l999) power spectrum is larger than our 
fiducial model by almost a factor of two, implying a sim- 
ilar correction factor to the global effect of minihaloes. In 
fact, we have conducted a comparison simu lation for the 'ex- 
treme suppression' case in which we use the lEisenstein fc Hul 
( 1999) power spectrum for the minihalo component, and we 
obtain an ionization fraction which is ~ 90% of our fiducial 
run. Since currently there are no observational constraints 
to distinguish one of these transfer functions from the other, 
we have chosen to remain consistent with our large-scale N- 
body simulations, but again, this could result in a conser- 
vative estimate of the minihalo photon consumption during 
reionization. 

A similar uncertainty arises from our choice of an ana- 
lytic form for the MH mass function. Over the relevant red- 
shifts , shifting from the Press-Schechter mass function to a 
ISheth fc Tormenl J2002I) form would decrease the total mass 



in minihalos by roughly a factor of 1.5. However, there is 
little evidence to motivate the use of this more complicated 
expression at very high redshifts. In fact, the few simulations 
of high-z small-scale structure formation that exist seem 
to agree reasonably well with the Press-Schechter approach 
JShapirdlioOll : Ijang-Condell fc Hernauisill200lt ICen et alJ 
12004) . 

The resolution of the simulations is not able to capture 
the IGM density distribution on very small scales. Usually 
this is taken into account, both in numerical simulations and 
semi- analytic approaches, through the so called gas clump- 
ing factor. As the main goal of this study is to discuss the ef- 
fect of MHs on the reionization process through a direct com- 
parison with the simulations described in CSW and CFW, 
we have not included this effect. Its inclusion would delay the 
reionization process of an amount depending on the adopted 
value or expression for the clumping factor (see e.g. ISS). It 
should be noted, however, that according to the results of 
ISS, any increase of the global photon consumption due to 
minihaloes is combined with the increased consumption due 
to recombinations in the clumpy IGM, often amplifying their 
impact. Any complete treatment of the effect of small-scale 
structures reionization should therefore include both effects. 



6 CONCLUSIONS 

We have investigated the impact of minihaloes on the reion- 
ization process by means of a combination of high-resolution 
N-body simulations (to describe the dark matter and diffuse 
gas evolution), high-resolution hydrodynamical simulations 
(to examine the photoevaporation of minihaloes), a semi- 
analytic model of minihalo formation (to follow their biased 
distribution), and the Monte Carlo radiative code CRASH 
(to follow the propagation of ionizing photons). We have 
studied the process assuming different parameters that regu- 
late the ionizing photon emission and different prescriptions 
for the evolution of minihaloes. The main results discussed 
in this paper can be summarized as follows. 

• If minihalo formation in a cell is completely sup- 
pressed once the first photon packet arrives at the cell, 
i.e. in partially or fully ionized cells ('extreme suppression' 
case) , we find that their effect on cosmic reionization is mod- 
est, with a volume averaged ionization fraction that is only 

15% lower than the one when minihaloes are ignored. 
Only less than 20% of the total emitted ionizing photons 
are absorbed by minihaloes at any redshift. Complete reion- 
ization is not delayed significantly by the presence of mini- 
haloes. 

• If minihalo formation is not suppressed ('no suppres- 
sion' case), then up to 80-90% of the emitted ionizing pho- 
tons are absorbed by MHs and complete overlap can be de- 
layed by as much as A z ~ 2. 

The 'no suppression' case is in good agreement with 
the results of the semi-analytical approach described in 
ISS. However, we fin d that the impact of MHs is sm aller 
than the estimat e s bv l Haiman. Abel fc Madau I i200(J) and 
Bark ana fc Loebl ([2002) , as these authors have overesti- 
mated the num ber of photons required to photoevaporate 
a MH. In fact, iBarkana fc Loebl 1120021) use static models 
of clouds in thermal and ionization equilibrium without 
accounting for gas dynamics. For this reason they fail to 
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capture some essential physics a nd overestimate the num- 
ber of recombinations inside MHs. lHaiman. Abel fc Madau I 
(2000) employ hydrodynamic simulations to study photo- 
evaporation, but they do n ot inc lude radiative transfer. 
IShapiro. Iliev fc Raeal J2004h and llliev. Shapiro, fc Raea I 
l|2005h (on which our calculations are based) found that the 
effect of ignoring radiative transfer and its feedback on the 
gas dynamics is to significantly overestimate the number of 
ionizing photons required to evaporate a MH. 

As discussed in § 15.31 there are still uncertainties in our 
model, some of which are conservative in terms of estimat- 
ing the photon consumption per minihalo. However, the 'no 
suppression' case certainly provides a liberal overestimate of 
the MH number density, which is likely to more than balance 
these effects. On the other hand, the 'extreme suppression' 
case provides a lower limit to the average MH number den- 
sity, since it assumes that no new MHs form in any cell 
that has ever been exposed to ionizing photons, no matter 
how weakly or long ago, and thus it is likely to bracket the 
MH correction effect from below. It is with some confidence, 
then, that we can limit the true impact of minihaloes to lie 
between the extreme cases considered here, although future 
investigations will be necessary to limit this range further. 
Thus, our results indicate that the photon sink provided 
by these structures is not sufficient by itself to force a long- 
lasting phase of partial IGM ionization. This agrees with the 
conclusion of ISS for the cases that neglected the small-scale 
clumping of the IGM. Thus the resolution of the apparent 
conflict between the high electron scattering optical depth in 
the WMAP data and the low reionization redshift inferred 
by QSO absorption line experiments should rely on addi- 
tional processes, such as the evolving small-scale IGM gas 
clumping factor, feedback effects or a transition in the star 
formation mode. 
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APPENDIX A: 
CELLS 



MINIHALOES IN SOURCE 



To account for the presence of MHs in computational cells 
that contain ionizing sources, we must take a slightly differ- 
ent approach from that described in §4. The reason is that 
at high redshift the ionizing sources were rare peaks of the 
density distribution, and thus minihaloes were strongly clus- 
tered around them. To account better for this clustering, we 
calculate the total number of photons absorbed by the gas 
(both diffuse and MHs) in that cell, N s , using a simplified 
version of the analytical model we developed previously in 
ISS. In the current simulations we group all sources within 
each cell into a single source with a mass equal to the sum 
of the individual masses. Note however that such multiple 
sources only occur at the latest redshifts simulated and thus 
this approximation has only a minor effect on our results. 

From ISS, the rate of expansion of the comoving volume 
about an ionized source of mass M 3 is given by 



dVi 
dt 



NM)- 1 - as (1 + zf n° H {Vi - V ) 



(Al) 

1- fcoll{Ms,r)+^rc(Ms,Z,F ) ' 

where n° H is the average comoving number density of hydro- 
gen, a_s = 2.6 x 10~ 13 cm 3 s _1 is the case B recombination 
coefficient for hydrogen at 10 4 K, Vb is the comoving volume 
initially carved out by the material making up the sources, 
and 

rAI, 

/coll, MH, arc (M s , Z, f) = 



, , dn(M,z, r\M 3 ) , , , 
dM — i — — — M , (A2) 



dM 



and 

£src(M s ,z,F ,r) 



dM x 



(A3) 



are the MH collapse fraction and average number of ex- 
tra absorbed photons by MHs at a Lagrangian distance r 
from a source of mass M 3 at a redshift z shining with a 
flux _Fb. Both f C oii{M s , z, r) and £ BTC (M S , z, Fo,r) make use 
of the biased number density of minihaloes, forming at a 
distance r from the source halo. To calculate this quantity, 
we make use of the analytical form alism des cribed in detail 
in IScannapieco fc Barkanal J2002T| (see also iPorciani et alJ 
1998j). Here the biased number density is given as 



(M,z,r\M 3 



dn 
dM 

where dn/dM : 



dn 1 



-(M,z,M s ,z,r) 



dn 



dMdM 3 y ' ' ' ' ' VdM s 
the usual PS mass function and 
dMdM {M,z,M s ,z,r) is the bivariate mass function that 
gives the product of the differential number densities at 
two points separated by an initial comoving distance r, at 
any two masses and redshift. Note that this expression in- 
terpolates smoothl y between all s tanda rd an alytical limits, 
includ ing those of iMo fc Whitel (ll996Tl and lLacev fc CoTel 
l|l993l) . a nd it has also been carefully v alidated against sim- 
ulations iScannapieco fc Thackerll2005Tl . 

Finally, we compute N s by integrating eq. (IAH up un- 
til the time t CG n at which the total mass contained within 
the ionized region is equal to the mass contained within the 
source cell. The number of photons absorbed is then simply 
N s = tN~f, which is imposed by adopting a lower effective 
escape fraction, such that the total number of ionizing pho- 
tons leaving the source cell is decreased. 



— (M,z,r\M 3 )MZ(M,z,F ), 



